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Abstract — We show that for continuous time dynamical sys- 
tems described by polynomial differential equations of modest 
degree (typically equal to three), the following decision problems 
which arise in numerous areas of systems and control theory 
cannot have a polynomial time (or even pseudo-polynomial 
time) algorithm unless P=NP: local attractivity of an equi- 
librium point, stability of an equilibrium point in the sense 
of Lyapunov, boundedness of trajectories, convergence of all 
trajectories in a ball to a given equilibrium point, existence of 
a quadratic Lyapunov function, invariance of a ball, invariance 
of a quartic semialgebraic set under linear dynamics, local 
collision avoidance, and existence of a stabilizing control law. 
We also extend our earlier NP-hardness proof of testing local 
asymptotic stability for polynomial vector fields to the case of 
trigonometric differential equations of degree four. 

I. Introduction 

Polynomial and trigonometric differential equations appear 
ubiquitously in a variety of application domains including 
robotics, economics, mathematical biology, and chemical 
engineering, among others [11]. The equations of motion 
for most robotic systems for example can be described by 
the familiar manipulator equations [19] which give rise 
to systems of differential equations that are a mixture of 
polynomial and trigonometric terms in the state variables. In 
mathematical biology and economics, polynomial differential 
equations such as the Lotka-Volterra model and its variants 
[14] are used to model population dynamics and competition 
among entities in an economy. The dynamics of many chem- 
ical processes [7] are also naturally modeled by polynomial 
differential equations. Aside from these specific examples, 
differential equations in numerous application domains are 
commonly approximated as polynomials. 

While mature computational tools exist for the numerical 
solution of such differential equations [5], in most of the 
application domains described above, one is not typically 
interested in particular solutions of the system. Rather, 
qualitative properties of the differential equations are of 
central importance. For example, one may be interested in 
the safety of a robot performing a certain dynamic task, or in 
determining if the population of certain species diminishes 
below a critical threshold. The former example is related 
to the stability of the control system employed by the robot 
while the latter can be addressed by defining an "acceptable" 
set of population numbers and asking whether this set is 
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invariant (i.e. if the populations start off in this set, will 
they always remain within the set?). In a similar vein, one 
can ask if trajectories in a model of epidemic spread remain 
bounded. 

The study of such qualitative properties of differential 
equations has an extensive literature and numerous algo- 
rithms have been proposed for addressing these questions 
computationally. However, for all but the simplest cases 
(e.g. the case of linear systems), the problems still lack 
satisfactory (i.e. exact and efficient) algorithms. This obser- 
vation motivates the study of the fundamental computational 
complexity of these problems in order to establish theoretical 
bounds on the efficiency of algorithms that attempt to answer 
these questions. Such complexity results may play an impor- 
tant role in shaping the search for practical algorithms for 
these problems by limiting the kinds of algorithms one can 
possibly hope to obtain. Further, by understanding exactly 
where the complexity of these problems stems from, we 
can hope to find approximations and relaxations that are 
more amenable to efficient solutions while still maintaining 
practical relevance. 

Questions of complexity related to qualitative properties of 
differential equations have long been of theoretical interest. 
A natural question one can ask is if the stability of a system 
of polynomial differential equations can be decided by a 
Turing machine in finite time. In [4], Arnold made a well- 
known conjecture that the contrary is true; i.e. the question is 
undecidable. To the authors' knowledge, even though some 
variants of the question have been studied and answered [12], 
[16], the question in its original form is so far unresolved. 
Although the results in this paper do not resolve Arnold's 
question, they provide lower bounds on the computational 
complexity of deciding local asymptotic stability and several 
similar and related problems. 

The primary challenge in establishing such complexity re- 
sults lies in relating the properties of the continuous solutions 
of polynomial and trigonometric differential equations to the 
combinatorial problems for which complexity results have 
been previously established. Explicitly mapping trajectories 
of a system (which typically one does not have access to 
exactly) to objects in combinatorial problems seems to be 
a hopeless approach. The main idea that allows us to by- 
pass this apparent challenge is to relate the combinatorial 
problem to properties of Lyapunov functions that prove 
stability /invariance of differential equations. All the results 
in this paper exploit this idea in one way or another 

The organization of the paper is as follows. After stating 
some preHminaries in Section |Il] we show in Section [Ul] 



that deciding local asymptotic stability for trigonometric 
polynomials of degree four is strongly NP-hard. While this 
result is an extension of the results presented in [2] (which 
proves the corresponding result for cubic polynomial vector 
fields), the decision problem is of independent interest, 
particularly in the field of robotics. This is due to the fact that 
most mechanical systems can be modeled by the manipulator 
equations, which result in vector fields whose degrees are 
dominated by trigonometric terms. In Section IIVI we prove 
that the following decision problems are strongly NP-hard 
for polynomial vector fields of degree d: 

• Invariance of a ball {d ~ 3), 

> Invariance of a basic semialgebraic set defined by a 
quartic polynomial {d = 1), 

• Inclusion of the unit ball in the region of attraction of 
an equilibrium point (d ~ 3), 

> Local attractivity of an equilibrium point (d = 3), 

« Stability of an equilibrium point in the sense of Lya- 

punov {d = 4), 
« Boundedness of trajectories (d = 3), 

• Existence of a quadratic Lyapunov function {d = 3), 

• Local collision avoidance {d = 4), 

• Existence of a stabilizing controller {d = 3). 

These notions are all formally defined in Section |IV] The 
input to these problems is an ordered list of coefficients 
(expressed as rational numbers) defining the polynomial or 
trigonometric vector field. Establishing NP-hardness of these 
problems implies that unless P=NP, it is not possible to 
provide an algorithm that can have a running time bounded 
by a polynomial in the number of bits required to represent 
the input. Further, all the NP-hardness results in this paper 
are in the strong sense (as opposed to weakly NP-hard 
problems like KNAPSACK or SUBSET SUM). This implies 
that the problems remain NP-hard even when the bit length of 
the coefficients (i.e. the input) is 0{log{n)) (here, n is the 
dimension of the state space). Unless P=NP, even pseudo- 
polynomial time algorithms cannot exist for strongly NP- 
hard problems; see [15] for more details and definitions. In 
particular, our results suggest that none of the numerous 
recent techniques for systems analysis based on convex 
optimization (e.g. in terms of linear programs, linear matrix 
inequalities, or sum of squares programs) can be exact, 
unless the size of the formulated optimization problems are 
exponential in the input. 

We refer the reader interested in computational complexity 
in systems and control to the outstanding survey papers [10], 
[16], [21], [8], [9] and references therein. 

II. A FEW PRELIMINARIES ON FORMS 

Many of the results in this paper will make use of homoge- 
neous polynomials. A multivariate polynomial p : K" — >■ R 
is homogeneous (of degree d) if it satisfies p{\x) = \'^p{x) 
for all X G M" and all A G R. This condition is equivalent to 
all monomials of p having the same degree. A homogeneous 
polynomial is also called a form. Observe that products 
of forms are again forms, and that the components of the 
gradient of a form are forms of one fewer degree. We will 



make frequent references to the following useful identity for 
homogeneous functions due to Euler: 

p{x) = —x^Vp{x). 

Here, p is a homogeneous function of degree d and Vp 
denotes its gradient vector. The identity is easily derived by 
differentiating both sides of the above equation with respect 
to A and setting A = 1. 

The degree of a polynomial vector field x ~ f{x), with 
/ : K" R", is defined to be the largest degree of 
the components of /. We say that the vector field / is 
homogeneous if all components of / are forms of the same 
degree. Finally, a form p : R" — 5- R is said to be positive 
definite if p{x) > for all nonzero x in R". 

III. Complexity of deciding local asymptotic 

STABILITY OF TRIGONOMETRIC VECTOR FIELDS 

In this section, we prove that deciding local asymptotic 
stability of trigonometric vector fields of degree four is 
strongly NP-hard. 

Definition 1: The zero equilibrium point of a dynamical 
system x = f{x) is stable in the sense of Lyapunov if V e > 
0, 3 (5 > such that ||2;(0)|| < S =^ \\x{t)\\ < e,\/t > 0. 
We say that the equlibrium point is locally asymptotically 
stable if it is stable in the sense of Lyapunov and there exists 
e > such that ||a;(0)|| < e limt^oo x{t) ~ 0. 

Theorem 3.1: Given a trigonometric vector field of degree 
four, it is strongly NP-hard to decide if it is locally asymp- 
totically stable. 

The proof will be via a reduction from ONE-IN-THREE 
3SAT, which is known to be NP-complete [15]. An instance 
of ONE-IN-THREE 3SAT consists of an expression made 
of conjunctions of clauses. Each clause is the logical OR 
of three literals, and each literal is either a variable or its 
negation. The problem is to decide if there is a boolean 
assignment of the variables that results in the expression 
being true and each clause having exactly one true literal. 

Given an instance of ONE-IN-THREE 3SAT in the vari- 
ables 6i, . . . , 6„, we first construct a degree 4 trigonometric 
polynomial. To avoid introducing unnecessary notation, we 
present this construction on a single instance of ONE- 
IN-THREE 3SAT. This example should elucidate how the 
procedure works in the general case. Given an instance of 
ONE-IN-THREE 3SAT such as: 

(6i V62V63)A(6i V64V65) 

we construct the following quartic trigonometric polynomial: 

t{x) = S,7=i sin2(a:;,)(l - sH^i)? 

+ [sin(a;i) + sin(a;2) + sin(a:;3) — 1]^ 

+ [sin(a;i) + (1 - sin(a;4)) + sin(x5) - 1]^ 

The first term in t{x) is not specific to the particular 
instance under consideration and will always appear in our 

'An assignment to a variable will denote an element from {0, 1} with 
coiTesponding to false and 1 coiTesponding to true. 



construction. The following terms are constructed by taking 
each clause appearing in the boolean formula and replacing 
a variable bi with sin(.Ti) and substituting + in place of V. If 
the negation of a variable appears in the clause, we replace 
it by 1 — sin(a:;i). Each resulting expression (corresponding 
to an individual clause) is subtracted by 1 and then squared. 

We then introduce a single new variable y and define z :— 
[x,y]'^ . We construct a new trigonometric polynomial: 

thiz) = S"^i sin^(a;i)(sin(2/) - sin(a;i))^ 

+ [sin(a:i) sin(j/) + sin(a;2) sin(y) . . . 

+ sin(a;3) sm{y) - sin^(2/)]^ 

+ [sin(xi) sin(j/) + {sm^{y) — sin(x4) sin(y)) . . . 

+ sin(a::5) sm{y) - sm^{y)f 

Here, we first replaced sin(a;i) with ^s'i^^(^y) ^'^'^ ^^^'^ multi- 
plied the entire resulting expression by sin^(?/). Observe that 
by doing so, th becomes a homogeneous function of sin(z). 

Definition 2: A function g(z) is locally positive definite if 
there exists e > such that ||z|| < e,z ^ ==^ g{z) > 0. 

Lemma 3.2: The trigonometric polynomial th{z) is lo- 
cally positive definite if and only if the ONE-IN-THREE 
3SAT instance it was derived from is unsatisfiable. 

Proof: To see this, first note that by construction 
t{x) and th{z) are sums of squares and hence nonnegative. 
Also, by construction, the only zeros of t{x) occur when 
sin(xj) e {0,1}. Now, suppose that the ONE-IN-THREE 
3SAT instance has a satisfying assignment of variables b*. 
Then, we observe that t{x) must have a zero if we substitute 
Xi = sin^^(&,*). Thus, by construction, t/i(z) has a zero 
when 

Xi = sin (6*) 
y = sin"^(l). 

Further, since is homogeneous in sin(z), for all a > 0, 
we have th{x{a), y{a)) = 0, where 

Xi{a) = sin~"'^(a6*) 
y{a) = sin~^(a). 

Note that as a ^ 0, {x{a),y{a)) gets arbitrarily close to 
the origin. Hence tji{z) is not locally positive. 

Now, to prove the converse, suppose the ONE-IN-THREE 
3SAT instance is unsatisfiable. We know that at least one 
of the clauses in the ONE-IN-THREE 3SAT instance must 
have either no true literals or more than one true hteral for 
all possible assignments to the variables. Thus, at least one 
of the terms of t{x) must always be positive. If t{x) does 
not have a zero, by construction, can only have a zero 
when sin(y) = 0. However, substituting sin(2/) = in tfi{z) 
results in the expression siii^(a;i). Thus, th{z) has a zero 
only when sin(x-i) ~ and sin(2/) ~ 0. Hence, in some 
neighborhood around the origin (— vr < < tt and — vr < 
y < tt), these equations cannot be satisfied and thus th{z) is 
locally positive. ■ 

Corollary 3.1: Checking local positivity of quartic 
trigonometric functions is strongly NP-hardo 

^Our reduction from ONE-IN-THREE 3SAT to checking local positivity 
is clearly polynomial in length. 



Le}mna 3.3: If the ONE-IN-THREE 3SAT instance is 
unsatisfiable, there exists a neighborhood around the origin 
in which \7th{z) does not vanish (except at the origin). 

Proof: Let s := [sin(.Ti), sin(.T2), . . . sin(j/)]^ and 
denote the function ts{s) to be the function th{z) viewed 
as a function of s. Note that ts is a quartic form in s. The 
expression for the gradient of t/i(z) with respect to z can be 
written as 

Vt/i(z) = diag(cos(z^))Vsts(s) 

where diag(cos(zi)) is a diagonal matrix with its diagonal 
elements set to cos(zi). Using Euler's identity for homoge- 
neous functions, we have ts{s) — js'^V sts{s). Thus, we 
see that when ts(s) is nonzero, Vsts{s) cannot vanish at 
that point. Since for a small enough neighborhood around 
the origin (—^<Xi<^ and — ^ < j/ < ^) we know that 
ts(s) is positive and cos(zi) is nonzero for all i, we have that 
Vth{z) also does not vanish in this neighborhood (except at 
the origin). ■ 

Corollary 3.2: Given a quartic trigonometric polynomial, 
it is NP-hard to decide if there exists a neighborhood around 
the origin where the function is positive and its gradient does 
not vanish (except at the origin). 

We now present a polynomial time reduction from the 
decision problem stated in Corollary 13.21 to the problem of 
checking local asymptotic stability of quartic trigonometric 
vector fields. This completes the proof of Theorem 13. II 

Proof: [Proof of Theorem 13.11 Given a quartic trigono- 
metric form th{z), we construct the following continuous 
time dynamical system: 

i = -~Vth{z). 

Observe that this a quartic trigonometric vector field. We 
claim that this system is locally asymptotically stable if and 
only if there is a neighborhood around the origin such that 
th{z) is positive definite and \7th{z) does not vanish except 
at the origin. To prove the claim, we start by noting that by 
construction tji{z) is always negative semidefinite: 

ihiz) = (Vi„(z),-Vi„(z)) = -\\Vt„iz)\\\ 

Suppose first that there is a neighborhood in which tji{z) 
is positive definite and Vtfi{z) does not vanish except at 
the origin. The condition on the gradient implies from 
the equation above that ih{z) is locally negative definite. 
This, together with local positivity of th, implies that th is 
a locally valid Lyapunov function for the system. Hence, 
local asymptotic stability follows from Lyapunov's stability 
theorem (see e.g. [18, p. 124]). 

To see the converse, suppose that the vector field is 
locally asymptotically stable. There is therefore an open 
neighborhood Bs around the origin, where the trajectories 
converge to the origin. We first observe that Vth{z) does not 
vanish in this neighborhood since if it did, z would equal zero 
at that point. This would contradict local asymptotic stability 
(since this point would be an equilibrium point and would 
not asymptotically converge to the origin). 



Next, we prove that th{z) is positive definite in Bs- 
Suppose first that there exists a point zq in Bs such that 
th{zo) < 0. Consider the trajectory that starts off at zq. 
Since ih{z) = — ||Vt/i(z)|p is non-increasing everywhere, 
we see that at all points in time, the value of th{z) remains 
negative and thus the trajectory cannot go to the origin (since 
th{0) ~ 0). This contradicts local asymptotic stability. To 
prove that t/i(z) is in fact strictly positive in Bs, suppose 
for the sake of contradiction that there is a point z* E Bs 
such that th{z*) = 0. Since tfi{z) was just shown to be 
nonnegative in Bs and Bs is an open set, z* is a local 
minimum. Thus, Vth{z*) = and z* is a fixed point. This 
again contradicts local asymptotic stability. ■ 

IV. Complexity of several qualitative properties 

OF POLYNOMIAL VECTOR FIELDS 

As we remarked earlier, NP-hardness of testing local and 
global asymptotic stability of polynomial vector fields of 
degree 3 has already been established in our earlier work [2], 
[1, Chap. 4]. In this section, we prove that deciding several 
other important properties of polynomial vector fields is also 
NP-hard. For many of these properties, our proof of NP- 
hardness builds on the proof in [2]. Whenever a property 
has to do with an equilibrium point, we take this equilibrium 
point to be at the origin. In what follows, the norm ||.|| is 
always the Euclidean norm, and the notation Br denotes the 
ball of radius r; i.e., B,. :~ | ||a;|| < r}. 

Theorem 4.1: For polynomial differential equations of de- 
gree d (with d specified below), the following properties are 
NP-hard to decide: 

(a) 0? = 3, Invariance of a ball: Va;(0) with ||a;(0)|| < 1, 

\\x{t)\\ < 1, Vt > 0. 

(b) d — 1, Invariance of a basic semialgebraic set defined 
by a quartic polynomial: \fx{0) £ S, 

x{t) eS, Vt> 0, 

where S := {x \ p{x) < 1} and p is a given form of degree 
four. 

(c) d = 3, Inclusion of a ball in the region of attraction of 
an equilibrium point: \/x{Qi) with ||a;(0)|| < 1, 

x{t) — > 0, as t ^ oo. 

(d) d = 3, Local attractivity of an equilibrium point: 36 > 
such that Vx(0) G Bs, 

x{t) — > 0, as t ^ CO. 

(e) d = 4, Stability of an equilibriym point in the sense of 
Lyapunov: Ve > 0, 3(5 = S{e) such that 

|a;(0)|| <S ^ \\x{t)\\ < e, > 0. 

(f) d = 3, Boundedness of trajectories: Va;(0), 3r — r{x{0)) 
such that 

\\x(t)\\ <r, yt> 0. 



(g) d = 3, Existence of a local quadratic Lyapunov function: 
35 > and a quadratic Lyapunov function V{x) = x^Px 
such that V{x) > for all x E Bs, x (or equivalently 
P ^ 0), and 

Vix) < 0, Va; e Bs, x ^ 0. 

(h) d ~ 4, Local collision avoidance: 35 > such that 

Vx(0) e Bs, 

x{t) ^S, Vt> 0, 
where 5 is a given poly tope. 

(i) d = 3, Existence of a stabilizing controller: There exists 
a particular (e.g. smooth, or polynomial of fixed degree) 
control law u{x) that makes the origin of 

X = f{x) + g{x)u{x) 

locally asymptotically stable, where / and g ^ Q here have 
degrees 3. 

Proof: (a) The proofs of parts (a) and (b) of the theorem 
are based on a reduction from the polynomial nonnegativity 
problem: given a (homogeneous) polynomial p : R" — > M, 
decide whether p(.t) > 0, Vx G R"? If the degree of p is even 
and larger or equal than 4, the problem is well-known to be 
NP-hard. This follows e.g. as an immediate consequence of 
NP-hardness testing matrix copositivity [20]. (The original 
proof of NP-hardness of checking matrix copositivity in [20] 
is via a reduction from the subset sum problem and only 
establishes weak NP-hardness. However, reductions from the 
stable set problem to matrix copositivity are also known [13] 
and they result in NP-hardness in the strong sense.) 

We now proceed with the proof of (a). Given a quartic 
form p, we construct the vector field 

X ~ -~Vp{x). 

Note that the vector field has degree 3 and is homogeneous. 
We claim that the unit ball Bi is invariant under the 
trajectories of this system if and only if p is nonnegative. 
This of course establishes the desired NP-hardness result. 
To prove the claim, consider the function V{x) := 
Clearly, Bi is invariant under the trajectories of i = —Vp{x) 
if and only if V{x) < for all x with ||a;|| = 1. Since V is 
a homogeneous polynomial, this condition is equivalent to 
having V nonpositive for all x G M". However, from Euler's 
identity we have 

V{x) = {2x, -Vp{x)) = ~8p{x). 

(b) Once again, we provide a reduction from the problem 
of checking nonnegativity of quartic forms. Given a quartic 
form p, we let the set S be defined as 5 = {x | p{x) < 1}. 
Let us consider the linear dynamical system 

X = —X. 

We claim that S is invariant under the trajectories of this 
linear system if and only if p is nonnegative. To see this, 
consider the derivative p of p along the trajectories of 



X ~ — X and note its homogeneity. With the same reasoning 
as in the proof of part (a), p{x) < for all x G M" if and 
only if the set S is invariant. From Euler's identity, we have 

p{x) = (Vp(x),— x) = — 4p(a;). 

Note that the role of the dynamics and the gradient of the 
"Lyapunov function" are swapped in the proofs of (a) and 
(b). 

A precursor to the proofs of (c)-(h). The proofs of (c)- 
(h) are all via reductions from the problem of testing local 
asymptotic stability (las) of cubic vector fields. This problem 
has been shown to be NP-hard in [2]. All of the reductions 
that follow are rather simple but subtly rely on the specific 
construction in [2]. For this reason, we need to first recall 
some elements of that construction. 

In a manner similar to that of the proof of Theorem l3.1l the 
reduction in [2] takes an instance of ONE-IN-THREE 3SAT 
and constructs a homogeneous polynomial V : M" — > M of 
degree 4 that is positive definite if and only if the ONE- 
IN-THREE 3SAT instance is not satisfiable. Then, a cubic 
vector field is constructed as follows: 

X -.^ fix) - -VV{x). (1) 

We recall the following facts from [2] about this vector field: 

1) The vector field in ([T]i is homogeneous and therefore 
it is locally asymptotically stable if and only if it is 
globally asymptotically stable. 

2) The vector field is (locally or globally) asymptotically 
stable if and only if V is positive definite. Hence, 
the system, if asymptotically stable, by construction 
always admits a quartic Lyapunov function. 

3) If the vector field is not asymptotically stable, then 
there always exists a nonzero point x G {0, 1}" such 
that f{x) = 0; i.e., a; is a nonzero equilibrium point. 

We now proceed with the proofs of (c)-(h). In what follows, 
f{x) will always refer to the vector field in (HJ. 

(c) The claim is an obvious implication of the homogeneity 
of /. Since f{Xx) = X^fix), for all A e M and all x G M", 
the origin is las if and only if for any r, all trajectories in 
Br converge to the originjj 

(d) If / is las, then of course it is by definition locally 
attractive. On the other hand, if / is not las, then f{x) = 
for some nonzero x G {0, 1}". By homogeneity of /, this 
implies that f{ax) = 0, Va > 0. Therefore, arbitrarily close 
to the origin we have stationary points and hence the origin 
cannot be locally attractive. 

(e) Let X* = (xf, . . . , xf^)'^ . Consider the vector field 

x = fix) + x\ (2) 

We claim that the origin of ^ is stable in the sense of 
Lyapunov if and only if the origin of ([T]i is las. Suppose 

^^For a general cubic vector field, validity of property (c) for a particular 
value of r is of course not necessary for local asymptotic stability. The reader 
should keep in mind that the class of homogeneous cubic vector fields is a 
subset of the class of all cubic vector fields, and hence any hardness result 
for this class immediately implies the same hardness result for all cubic 
vector fields. 



first that ([T]) is not las. Then we must have f{ax) = 
for some nonzero x G {0, 1}" and Va > 0. Therefore for 
the system (O, trajectories starting from any nonzero point 
on the line connecting the origin to x shoot out to infinity 
while staying on the line. (This is because on this line, the 
dynamics are simply i = a;'*.) As a result, stability in the 
sense of Lyapunov does not hold as there exists an e > (in 
fact for any e > 0), for which $S > such that trajectories 
starting in Bs stay in B^. Indeed, as we argued, arbitrarily 
close to the origin we have points that shoot out to infinity. 

Let us now show the converse. If ([T]i is las, then V is 
indeed a strict Lyapunov function for it; i.e. it is positive 
definite and has a negative definite derivative — ||Vy(x)|p. 
Using the same Lyapunov function for the system in (|2|i, we 
have 

V{x) = -\\VVix)\\' + {VV{x),x'). 

Note that the first term in this expression is a homogeneous 
polynomial of degree 6 while the second term is a homo- 
geneous polynomial of degree 7. Negative definiteness of 
the lower order term implies that there exists a positive real 
number S such that V{x) < for all nonzero x G Bs- 
This together with positive definiteness of V implies via 
Lyapunov's theorem that (|2| is las and hence stable in the 
sense of Lypunov. 

(f) Consider the vector field 

x = f{x) + x. (3) 

We claim that the trajectories of (O are bounded if and only 
if the origin of ([Hi is las. Suppose (HJ is not las. Then, as in 
the previous proof, there exists a line connecting the origin to 
a point X G {0, 1}" such that trajectories on this line escape 
to infinity. (In this case, the dynamics on this line is governed 
by i = X.) Hence, not all trajectories can be bounded. For the 
converse, suppose that ([T]) is las. Then V (resp. — ||Vy(a:;)|p) 
must be positive (resp. negative) definite. Now if we consider 
system ([3]), the derivative of V along its trajectories is given 
by 

V{x) = -||Vy(x)|p + {VV{x),x). 

Since the first term in this expression has degree 6 and the 
second term degree 4, there exists an r such that < for 
all X ^ Br- This condition however implies boundedness of 
trajectories; see e.g. [17]. 

(g) If / is not las, then there cannot be any local Lyapunov 
functions, in particular not a quadratic one. If / is las, then 
we claim the quadratic function W{x) = \\x\\'^ is a valid 
(and in fact global) Lyapunov function for it. This can be 
seen from Euler's identity 

W{x) = (2a;, -VV{x)) = ~8V{x), 

and by noting that V must be positive definite. 

(h) We define our dynamics to be the one in (O, and the 
polytope S to be 

n 

S = {x \ Xi>0,l <^Xi < 2}. 

1=1 



Suppose first that / is not las. Then by the argument 
given in (e), the system in (O has trajectories that start out 
arbitrarily close to the origin (at points of the type ax for 
some X e {0, 1}" and for arbitrarily small a), which exit 
on a straight line to infinity. Note that by doing so, such 
trajectories must cross S; i.e. there exists a positive real 
number a such that 1 < X^ILi '^^i — 2- Hence, there is 
no neighborhood around the origin whose trajectories avoid 
S. 

For the converse, suppose / is las. Then, we have shown 
while proving (e) that (|2|i must also be las and hence stable 
in the sense of Lyapunov. Therefore, there exists S > such 
that trajectories starting from Bg do not leave Bi — a ball 
that is disjoint from S. 

(i) Let / be as in ([T]i and g{x) ~ [xix"^ — a;^.T2)ll"^, 
where 1 denotes the vector of all ones. The following 
simple argument establishes the desired NP-hardness result 
irrespective of the type of control law we may seek (e.g. 
linear control law, cubic control law, smooth control law, 
or anything else). If / is las, then of course there exists 
a stabilizing controller, namely u = 0. If / is not las, 
then it must have an equilibrium point at a nonzero point 
X G {0,1}". Note that by construction, g vanishes at all 
such points. Since g is homogeneous, it also vanishes at all 
points ax for any scalar a. Therefore, arbitrarily close to 
the origin, there are equilibrium points that the control law 
u{x) cannot possibly remove. Hence there is no controller 
that can make the origin las. ■ 

Remark 4.1: Arguments similar to the one presented in 
the proof of (i) above can be given to show NP-hardness 
of deciding existence of a controller that establishes several 
other properties, e.g., invariance of the unit ball, inclusion of 
the unit ball in the region of attraction, etc. In the statement 
of (h), the fact that the set 5 is a polytope is clearly arbitrary. 
This choice is only made because "obstacles" are most 
commonly modeled in the literature as polytopes. We also 
note that a related problem of interest here is that of deciding, 
given two polytopes, whether all trajectories starting in one 
avoid the other This question is the complement of the usual 
reachability question, for which claims of undecidability 
have already appeared [16]; see also [6]. 

V. Conclusions and Open problems 

Under the assumption P^NP, we have shown the impossi- 
bility of polynomial time (or even pseudo-polynomial time) 
algorithms for ten decision problems that ubiquitously arise 
in control theory and the study of continuous time dynamical 
systems. Although our hardness results are valid even for 
very restricted classes of systems (e.g. gradient systems), it 
is of course still possible that these decision problems admit 
polynomial time algorithms for other special (and possibly 
important) subclasses of polynomial or trigonometric vector 
fields. 

Aside from extending the results of Theorem 14.11 to other 
classes of vector fields (such as trigonometric ones), the 
obvious class of questions that our work leaves open is to 
investigate the decidability of the decision problems studied 



in Theorem 14. 11 or their NP-hardness for polynomial vector 
fields of degree one or two. Although for linear systems 
some of these questions become easy, we expect that our 
hardness results can be strengthened to the case when the 
degree is 2. Quadratic vector fields already demonstrate very 
complex behaviour; for example, their stability does not 
imply existence of a polynomial Lyapunov function of any 
degree [3]. In general, one can reduce the degree of any 
vector field to two by introducing polynomially many new 
variables (see [16]). However, this operation may or may not 
preserve the property of the vector field which is of interest. 
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